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SUMMARY 


This  report  presents  the  o^jithematlcal  derivation  of  two  methods  for 
computing  the  radioactive  fallout  at  a  point.  Information  is  given  which 
will  permit  the  final  calculation  of  dose  or  dose  rate. 
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I .  IMTROJXJCTION 

It  1b  imperative,  in  the  nuclear  age,  that  we  be  able  to  evaluate  the 
effects  of  radiological  fallout.  In  order  to  obtain  an  adequate  repreBentation 
of  the  fallout  field  \dilch  may  result  from  an  atomic  campaign,  it  is  necessary 
first  to  understand  the  mechanism  of  fallout;  second,  to  have  a  simple 
realistic  method  for  estimating  the  fallout  field  from  a  single  bomb;  and, 
third,  to  have  a  simple  and  realistic  method  for  combining  the  patterns  from 
many  bombs  into  a  conq>rehen6ive  picture  of  the  distribution  and  intensity  of 
fallout. 

(l  2) 

Ihe  basic  mechanism  of  fallout  has  been  understood  for  some  time.'  *  ' 
However,  the  development  of  a  simple  schene  for  estimating  the  fallout  ftom 
a  single  bomb  has  been  delayed  because  of  the  uncertainties  of  the  distri¬ 
bution  of  activity  within  the  cloud  formed  by  the  detonation,  the  distri¬ 
bution  of  activity  with  particle  size,  the  fall  rate  of  the  pEurtlcles,  and 
the  determination  of  the  true  velocity  of  the  wind.  A  recent  study,  reported 
in  a  series  of  papers, addressed  the  problems  of  estimating  the 
necessary  parameters  and  determining  the  effect  of  changes  in  these  parsp- 
meters  on  the  final  fallout  patterns.  The  purpose  of  this  paper  is  to 
present  a  brief  outline  of  the  fallout  mechanism,  to  derive  from  this 
mechanism  a  mathematical  model  to  represent  the  main  featxures  of  the 
fallout  process,  and  to  provide  a  detailed  exposition  of  the  simple 
solutions  to  this  mathematical  model. 

It  is  expected  that,  having  derived  these  solutions  of  mathematical 
models  vdilch  can  adequately  describe  fallout  patterns  over  a  fairly  wide 
range  of  yield  and  wind  conditions,  the  next  step  can  be  taken;  namely. 
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the  dereloinient  of  easy  method*  for  applying  the  solutions*  It  has  long 
been  the  ala  of  those  vorking  on  the  theory  of  fhllout  to  derelop  a 
relatively  streamlined  iiietVvv>«l  for  coiqinitlng  fallout,  cmd  it  is  believed  that 
the  development  described  here  vlll  lead  to  such  a  tool. 
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II .  BASIC  FACTORS 

llie  detonation  of  a  nuclear  device  on  or  above  the  ground  creates  a 
hot  bubble  of  gas  containing  all  of  the  fission  products  of  the  detonation, 
nils  hot  bubble  will  rise  and,  because  of  the  strong  temperature  gradient, 
form  a  vortex  ring.  The  vortex  ring,  by  virtue  of  Its  stability,  will 
trap  the  fission  fragments  and  carry  them  aloft  Into  the  mushroom  part  of 
the  cloud.  If  the  detonation  Is  on  or  very  close  to  the  surface  of  the 
earth,  large  amounts  of  dirt  will  be  swept  Into  the  clrculatlcc  and  the 
fission  fragments  will  be  dex>oslted  on  peortlcles  of  veirlous  sizes.  If 
the  detonation  Is  so  high  In  the  air  that  the  vortex  ring  Is  formed  before 
appreciable  quantities  of  earth  are  swept  Into  It,  the  particles  which  are 
formed  will  be  extremely  small  and  will  not  fall  during  the  first  few  hours. 
The  distribution  of  the  activity  with  peurticle  size  has  a  direct  beeurlng 
on  the  close-in  feQJ.out.  If  there  Is  considerable  earth  swept  Into  the 
cloud,  approximately  80  per  cent  of  the  fission  fragments,  and  also  any 
induced  radioactivity,  will  be  lodged  on  particles  which  will  fall  to  the 
ground  within  24  hours. 
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III .  IHE  MODELS 


This  paper  will  deal  only  with  weapoDB  In  the  megaton  range  which  are 
detonated  on  the  surface  of  the  earth.  This  choice  of  problem  limitation 
1b  made  because  it  represents  the  most  interesting  case  for  close-in 
fallout.  Smaller  yields  or  bxirsts  at  higher  elevations  will  produce 
smaller  and  less  Intense  patterns. 

The  conditions  of  the  cloud  from  5  to  7  min  after  the  burst  provide 
the  starting  point  for  the  fallout  calculations.  At  this  time  the  cloud  has 
stopped  rising,  the  vortex  ring  circulation  has  slowed  considerably,  and  the 
cloud  has  cool§d.  Ibe  fission  fl-agments,  induced  activity,  and  earth  have 
formed  radioactive  particles  which  are  distributed  throughout  the  cloud. 

For  purposes  of  csilculatlon,  the  cloud  can  be  described  by  the,  distribution 
of  actlvily  as  a  function  of  position  in  the  cloud  and  the  size  of  particles 
on  ^ich  the  activity  is  lodged.  The  subsequent  history  of  these  particles 
is  determined  by  the  fall  velocity  and  the  wind. 

If  the  winds  are  known  as  a  function  of  time  and  spac'  ,  and  the  rate 
of  fedl  of  the  particles  can  be  determined,  then  the  final  ground  position 
of  a  particle  of  radius  r,  and  initial  height,  H,  can  be  determined.  This 
can  be  expressed  mathematically.  If  1D  (r,H)  is  the  horizontal  distance 
^Ich  the  particle  travels  and  V  (x,y,h,t)  is  the  wind  velocity, 

T 

Id  (r,H)  =  /V  (x,y,h,t)  dt 
o 

where  the  Integration  follows  the  path  of  the  particle.  The  vertical 
velocity  of  the  particle  is  given  by 


W(r,H) 


dh 

3t 
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BO  that  the  horizontal  distance  traveled  is 

)D  (r,H)  -  I  I  dh  (1) 

Even  if  it  were  possible  to  determine  the  exact  winds  and  fall  velocities, 
the  complete  mathematical  expression  would  make  Eq^>  (l),  which  we  will  call 
the  transport  equation,  iii^>osslble  to  solve,  ^e  transport  eqiiatlon  will 
therefore  be  approximated  by  simple  functions  \dilch  are  in  accord  with  our 
state  of  knowledge  of  the  wind  and  fall  velocities. 

nie  "transport  eq.uatlon"  can,  in  principle,  be  Inverted  to  determine 
the  initial  position  and  size  of  the  particles  which  will  land  at  a  given 
point.  If  the  activity  is  Integrated  over  all  such  particles,  the  radiation 
from  this  point  can  be  determined.  Since  the  initial  cloud  condition 
is  expressed  as  a  density  function,  and  since  the  radloactlvl'ty  is  not 
affected  by  the  fall  and  transport  of  the  particles,  it  is  possible  to 
deal  only  with  the  fraction  of  actlvily  and  time  of  arrival.  The  Inltleil 
distribution  of  activity  in  the  cloud  may  be  expressed  as  A(r,H,R),  where  r 
is  the  radius  of  the  particles,  H  is  the  height  above  the  grotind,  and  R  is  the 
radial  dimension  of  the  cloud.  If  circular  symmetry  is  assumed,  A(r,H,R) 
is  the  fraction  of  activity  per  micron-kft-sq  n  mi. 

Now,  the  radioactivity  which  will  affect  people  is  a  function  of  the  den¬ 
sity  of  materlsJL  on  the  surface;  therefore,  if  the  activity  function  is  inte¬ 
grated  over  the  peurtlcle  sizes,  from  all  heists,  which  land  at  a  point,  the 
result  will  be  the  density  function  per  sq  mi. 

F(x,y)  =  J*  J  A(r,H,R)  dr  dH  (2) 

Hi 

This  will  be  termed  the  "fraction  equation"  emd  is  the  second  necessary 
equation  for  the  mathematical  model  of  fallout. 
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In  order  to  make  these  two  equations  tractable  It  is  necessary  to  make 
some  simplifying  assumptions.  The  assumptions  will  be  listed  in  two  groiips, 
without  any  detailed  rationale.  The  first  group  of  assuii5)tlons  are  necessary 
to  provide  any  type  of  solution  and  the  second  group  cure  needed  to  provide 
practical  computing  procedures. 

The  first  assumption  is  that  the  particles  fall  with  the  terminal 
velocity  of  spheres.  Hie  experimental  work  reported  by  Rapp  and  Sartor^^^ 
is  the  principal  Justification  for  this  assumption.  The  next  assumption  is 
that  the  winds,  as  measured  by  normal  meteorological  sounding  methods, 
provide  sufficient  detail  to  estimate  the  displacement  of  the  falling 
particles.  Essentially,  this  assumes  that  turbulence  will  have  no  effect 
on  the  particles.  This  problem  is  discussed  by  Rapp.^^^^  With  these  two 
restrictions  the  transport  equation  has  been  solved  numerica4-.y  for  space- 
and  time-varying  winds. 

In  order  to  solve  the  fraction  equation  the  distribution  function, 
A(r,H,R),  is  needed.  This  was  assumed  to  be  the  product  of  three  Independent 
distributions 

A(r,H,R)  «  A^(r)  AgCn) 

.A2(H)  and  Aj(R)  were  determined  empirically  ffom  observations  of  the 
distribution  of  activity  in  a  cloud  at  7  mln.^^^  was  assumed  to  be  a 

log  normal  distribution  and  the  parameters  of  this  distribution  were  chosen 
to  fit  the  fallout  pattern  of  a  test  shot,  subject  to  the  assumption  that 
AgCH),  Aj(R),  V  (x,y,h,t),  and  w(r,H)  were  all  known. The  limits  of 
integration  were  then  determined  from  the  transport  equation.  This  system 
of  equations  was  solved  numerically  for  the  wind  data  from  a  nuclear  test 
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vhlch  vas  not  used  in  deriving  the  model  amd  the  results  were  found  to  be 
adequate* 

It  Is  not  claimed  that  the  model  presented  here  is  unique  or  precise* 
The  model  Is,  however,  based  on  accepted  physical  principles,  is  consistent 
with  meteorological  practice,  is  cosq^able  with  model  experiments,  and  is 
not  inconsistent  with  the  results  of  nuclear  test  measurements,  nie 
solution  to  the  problem  of  fallout  conq>utatlcn  based  on  the  restrictions 
outlined  above  was  possible  only  by  a  lengthy  numericsil  Integratioa  using 
high-speed  cooqputlng  machinery.  In  order  to  provide  a  quicker  and  more 
sl]iq>le  scheme  for  making  fallout  cooputations  some  further  assumptions 
are  needed. 

If  a  mean  value  theorem  is  applied  to  the  velocity  term  in  the 
tiransport  equation,  then 

H 

ID  (r,H)  -  V  f  -  V  T(r,H)  (5) 

where  T(r,H)  is  the  time  it  takes  for  a  particle  of  radius  r  to  fall 
ttam  a  height,  H.  By  numerical  integration  of  the  terminal  velocities  of 
particles  from  various  heights,  it  was  found  that  the  fall  times  could  be 
approximated  quite  well  by 

T(r,H)  -  a(r)  ♦  (h) 

r 

The  functions  a(r)  and  B(r)  are  peuraffleters  which  vary  only  slowly  over  the 
range  of  r  and  therefore  mjiy  be  assumed  consteuit  over  liad.ted  ranges  of  r 
without  committing  any  serious  error.  When  using  the  transport  equation 
in  the  form  of  Eq.  (3)  the  appropriate  velocity  should  be  some  sort  of  a 
weighted  mean  velocity. 
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The  fraction  equation  is  further  simplified  by  replacing  the 

continuous  empirical  distributions  A2(H)  and  Aj(R)  by  a  constant  value 

between  limits.  Hius,  if  A^(R)  is  considered  to  be  a  constant,  Aj,  for 

0<R<R  ,  idaere  R  is  an  effective  radius,  it  can  be  taken  outside  the 
e  e 

integral  sign.  Similarly  AgCH)  may  be  replaced  by  a  constant,  and  tsdsen 
outside  the  integral  sign.  The  product  of  suid  A^  is  simply  the  reciprocal 
voltone  of  the  cloud,  so  that  the  fraction  equation  becomes 

rgCH) 

F(x,y)  =  f  A  (r)  dr  dH  (5) 

The  geometry  of  the  sittiatlon  is  shown  in  Fig.  1.  The  shaded  area 
represents  a  cross  section  of  the  cloud  and  is  assumed  to  have  a  constemt 
concentration  of  activity  in  space.  The  particle  of  radius  r  at  height  H 
at  a  distance,  R,  for  the  center  will  fall  to  the  ground  during  time  T 
(E3q.  4)  and  will  travel  a  distance,  |^|T,  where  V  is  the  mean  wind.  The 
height  of  the  base,  H^,  the  height  of  the  top,  H^,  and  the  effective  radius, 
define  the  limits  of  the  cloud.  The  coordinate  axes  are  H,  vertically 
upward;  x,  from  ground  zero  in  the  direction  of  ^  and  y  from  ground  zero 
perpendicular  to  the  direction  of 

The  direction  and  magnitude  of  ^  can  be  computed  by  many  methods. 

One  of  the  most  convenient  methods  is  to  add  vectorially  T/i6  of  the  500  mb 
wind  vector  (l8,000  ft),  5/l6  of  the  100  mb  wind  vector  (50,000  ft),  and 
U/l6  of  the  50  mb  wind  vector  (68,000  ft). 

There  is  one  final  bit  of  Information  needed  to  complete  the  method 
of  calculation:  the  variation  of  the  cloud  parameters  with  yield.  The 
Summary  Report  of  RAIID  Work  on  the  AFSV/P  Fallout  Project 


presents 
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some  scaling  laws  which  were  based  on  the  study  of  many  observed  clouds. 
Ihese  laws  were  the  result  of  the  fitting  of  least  squares  to  the  observa¬ 
tions.  Since  the  fitting  procedure  resulted  in  irrational  numbers  for 
the  powers  of  the  yield  of  the  device,  a  consistent  and  sirqple  set  of 
scaling  laws  is  desired.  It  is  assumed  that  the  effective  volume  of  the 
cloud  is  directly  proportional  to  the  blast  yield 


Vol  -  C^W 

^ere  W  is  the  yield  of  the  device.  Simple  expressions  for  the  effective 
radius  and  height  difference  were  found  to  be 


“t  -  “b  ■ 


was  found  to  be  22  kft,  and  C2  was  determined  as  8.O6  n  mi,  so  that 
2 

4^00  n  ml  kft.  In  addition  to  the  scaling  of  the  cloud  dimensions 
the  cloud  height  was  found  to  scale  as 

Hp  =  28.8 

In  all  these  scaling  lavs  horizontal  dlmensloDS  are  measvired  In  n  ml, 
vertical  dimensions  in  kft,  and  yield  in  KF. 
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IV.  HEIGHOyf  ARTICLE  SIZE  MmOD  OF  IimCRATION 

nie  purpose  of  this  section  Is  to  find  a  solution  to  Eq,.  (^),  '^Ich  will 
utilize  Eq.  (U),  to  determine  the  limits  of  Integration.  The  problem  now 
becomes  one  of  determining  the  range  of  particle  sizes  \dilch  will  reach  each 
point  downwind  of  ground  zero.  At  each  height,  H,  within  the  cloud,  the 
maximum  sized  particle,  r2^  which  lands  at  a  given  point  will  fall  from  the 
front  of  the  cloud  and  the  minimum  sized,  r^,  from  the  rear  of  the  cloud. 

If  the  cloud  were  infinitely  thin  the  activity  at  the  point  would  be  that 
fraction  of  the  activity  on  the  particles  in  the  range  from  r^^  to  Tg. 

However,  because  the  cloud  Is  of  finite  thickness,  the  range  of  particles 
Is  a  function  of  height  and  the  actlvl'ty  at  the  point  becomes  the  sum  of 
the  flractions  of  the  activity  on  the  particles  in  the  range  from  r^^  to  rg 
from  each  height  within  the  cloud. 

In  order  to  Integrate  Bq.  (3),  the  fonctlons  r^(H)  and  r2(H)  must  be 
determined.  The  radius  of  a  particle,  r,  falling  at  an  arbitrary  point 
on  the  ground  (x,y),  trom  an  Initial  height,  H,  throu^  a  constant  wind 
field,  ^  ,  is  uniquely  determined  by  Its  time  of  fall,  T(r,H).  It  has 
been  assumed  that  the  time  of  fall  may  be  approximated  by  Eq.  (4).  Let 
r^  and  rg^  be  the  radii  of  the  smallest  and  largest  particles  falling  from 
the  base  of  the  cloud  to  the  point  (x,y).  The  particles,  r^(H)  and  rg(H), 
which  are  vertically  above  and  rgg,  will  reach  the  same  point  or.  the 
groimd  at  the  same  time.  Thus 


RH  ^  B 

T(r,H)  =  a  +  ^  =  a  +  — ^ 

(6a) 

(6b) 
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Similarly 

ri(H)  - 


(6c) 


Here  it  has  been  assumed  that  a  and  S  remain  constant  over  the  range  of 
particles  from  the  bottom  to  the  top  of  the  cloud. 

The  activity-particle  size  distribution  has  been  approximated  by  a  log¬ 
normal  distribution  vlth  mean,  \i,  and  standard  deviation,  o  (u  ■  3>8  and 
o  ■  0.69)^  '  Using  the  log-normal  distribution  and  Eq.  (6)  for  the  limits 
of  Integration,  Eq.  (9)  becomes 


F(x,y)-^J  _  e 

Inr^^/S/m 


d(lnr)dH  (7) 


If  8  «  ,  and  ds  »  ^  d(lnr) 

then  F(x,y)  -  /  /  e'^/^  8® 

""b  «1 


dsdH 


vhere 


In  (r^VH^)  -  4 


■  X,  /  {  g  ,e, 

■  -  ^i} 

vhere  0  Is  the  standard  normal  distribution. 
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By  conqpletlng  the  square  In  the  exponent 

2(u  -  In  r^//C)  r  2ot 

■  2ce  -25- 

-  2oe^^^  “  ^  f  dt 

"b 

If  w  >  t  —  2o,  then  dt  ■  dv 
and 


An  Inspection  of  Pj^  shows  that  as  — >  00  Pj^  — H 

Pj/dH  — ^  1  as  — >  00 .  Uie  function  F^/£E  was  computed  for  several 

values  of  r^  and  plotted  on  log-probability  paper.  Ihe  results  indicated 
that 


j,  -  Hp.  Iherefore 


(10) 


With  the  choice  of  scaling  laws  discussed  in  Section  III,  the  function 

P. 

^  (r^)  iB  independent  of  yield,  ^(r^^)  is  then  a  log-nonnal  distribution 

with  mean  equal  to  3.63  and  0^  >  O.69. 
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This  new  log-normal  distrlhution  has.  In  effect,  reduced  the  problem 
to  one  of  cojiputing  fallout  from  an  infinitely  thin  cloud.  The  only  change 
in  the  distribution  is  foiind  in  the  mean,  vAiich  has  decreased.  This 
decrease  of  the  mean  represents  an  increase  of  activity  on  the  particles 
falling  from  the  base  by  an  amoimt  equal  to  the  sum  of  the  activity  on  the 
particles  falling  from  all  heights,  thus,  in  effect,  concentrating  all  the 
activity  at  the  base.  The  equation  for  fraction  down  becomes 


F(x,y)  =  (r2B)  )]•  (ll) 

TTR  L  J 

e 

Up  to  this  point  no  mention  has  been  made  of  the  meteorological 
variable,  wind  velocity,  with  the  exception  that  it  is  connected  to  the 
fraction  down  throiigh  the  fall  velocity  of  the  particles.  As  was  mentioned 
earlier  the  radius  of  each  particle  falling  from  height  H  determines  the 
time  down.  ®ils  time  determines  the  horizontal  distance  the  particle  will 
travel  in  a  given  wind  velocity.  In  order  to  solve  Eq.  (ll)  one  must  first 
find  rgg  and  r^. 

Figure  2  shows  a  cross  section  of  the  cloud  in  the  direction  of  the 


mean  wind.  This  figure  clearly  shows  that  the  particle  r^.  travels  a 
horizontal  distance,  x  —  R  ,  and  the  particle,  r,-,  a  horizontal  distance, 

X  +  R^.  With  a  mean  wind  of  ^  the  times  of  fall  of  these  particles  would  be 


X  -  R 

Tp  =  - ^  (12) 

V 

and 

X  +  R 

T,  =  -  _  --  (13) 

V 


This  gives  the  times  of  fall  for  the  largest  and  smallest  particles  on  the 
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x-axls.  Figure  5  shows  the  cloud  ftrom  the  top.  It  can  bo  seen  that 
particles  falling  from  the  cloud  at  a  distance,  y,  from  the  x~€ucls  will 
travel  different  horizontal  distances.  Ibese  distances  depend  upon  the 
effective  radius  of  the  cloud  at  y.  By  geometry  the  effective  radius 
becomes 

Ihus,  for  off-axis  points,  Eqs.  (12)  and  (13)  become 


and 


X  -  X  -  Rgf(y) 

■  ■I  I  I  »  a  ■  mi ■  — 

V  ? 


X  +  R  • 
_ e 

7 


X  +  Rgf(y) 


(12a) 


(15a) 


Here  f(y)  Is  Introduced  to  scale  the  cloud  radius  for  the  off-axis 
points.  By  geometry  f(y)  Is 


In  theory,  r^^^  and  con  be  computed  for  various  yields  from  Bi*  (4) 
and  the  scaling  lav  for  H^.  Because  of  the  dependence  of  a  and  p  on  r 
It  Is  convenient  to  graph  r  as  a  function  of  T  and  W  (see  Fig.  4). 

This  Biethod  may  also  be  used  to  ccmpute  the  fraction  down  at  points 
under  the  cloud  and  behind  ground  zero.  Ihe  range  of  particle  sizes 
falling  from  the  base  to  a  point  under  the  cloud  will  Include  all  particles 
larger  than  the  size  of  those  falling  from  the  rear  of  the  cloud.  Ihe 
function  ^  (^2B^  vlU  then  become  \inlty  and  the  fraction  down  equation 


becomes 


YwM  (MT) 
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Tim#  (hr) 


Fig.  4 — Particle  size  os  o  function  of  yield  and  time 
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F(x,y)=4-  (15) 

R 

e 

93ie  fraction  down  is  obviously  a  function  of  yield  and  wind  velocity. 

It  mis^t  be  helpful  at  this  point  to  describe  how  the  fraction  changes  with 
changing  yield  and  wind  velocity.  From  Eq.  (ll)  it  c£ui  easily  be  seen  that 
F(x,y)  varies  inversely  with  the  square  of  the  cloud  radius.  Thus,  as  the 
yield  Increases,  the  fraction  decreases. 

Another  change  in  the  fraction  can  be  attributed  to  the  cloud  radius. 

It  is  clear  from  Eqs.  (12)  and  (15)  that,  for  a  given  V,  the  period  during 
>dilch  particles  arrive  at  a  point  Increases  with  increasing  cloud  reuilus. 

The  resulting  Increase  in  the  range  of  particle  sizes  vlU  increase  the 
fraction  down.  In  connection  with  Eqs.  (12)  and  (l5)#  it  can  be  seen  that, 
for  a  given  (x,y),  an  Increase  of  wind  velocity  would  produce  a  decrease 
in  the  range  of  particle  sizes  and  a  shift  of  that  range  to  larger  particle 
sizes.  Hie  last  effect  is  attributed  to  the  increasing  Hg  with  increasing 
yield.  For  a  given  (x,y)  an  increase  in  yield  produces  a  shift  in  the 
range  of  particles  toward  the  large  end  of  the  spectrum.  All  these  effects 
working  together  make  any  method  of  scaling  the  fraction  for  yield  too 
cos^llcated  to  be  useful. 

In  the  next  few  paragraphs  a  step-by-step  procedure  for  computing 
the  fallout  at  a  point  will  be  outlined. 

The  quantities  given  for  fallout  patterns  are  yield  and  wind  velocity. 
As  an  example,  consider  a  U-l-lT  device  and  a  mean  wind  of  10  kn.  The 
procedure  is  as  follows  for  points  on  the  x  axis: 

A.  Determine  the  cloud  radius  and  l/ ttR^  from  Fig.  5 

R  =  15  n  mi 
e 

I/ttRq  =  O.OOlUl 
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Yield  (MT) 


Fig.  5  — Cloud  radius  (R,)  and  the  reciprocol  of  cloud  oreo^^^) 

os  0  function  of  yield  * 
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B.  Determine  the  dovnvlnd  distances  desired  and  record. 

C.  Compute  T2  and  from  Eqs.  (12)  and  (l3)  and  record. 

D.  Determine  the  particle  sizes  for  the  appropriate  yield  and 
times  from  Fig.  4  and  record. 

E.  Determine  ^  ^  record. 

P.  Take  the  difference  of  ^  (^23^  i'  (rgg),  then  multiply 

by  l/rrR^.  Hie  result  is  F(x,y). 

For  points  off  the  x-axls  the  same  procedure  is  foUoved  with  the  exception 
that  the  function  f(y)  is  used  in  determining  the  times  T2  and  T^. 

For  exanqiile,  if  x  «  30  and  y  *  l/SR,  then 

I,  .  .  3.7  hr 

10 

.  50  .  1?  Vl  -  (1/a)!  .  6.5  hr 


Using  these  times,  the  following  yalues  of  particle  size  are  obtained 
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Fig.6 
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V.  TIMB-^ARTICLE  SIZE  MEIHOD  OF  INTBGRATIOH 

We  proceed  to  perform  the  Integration  of  the  fraction  equation  to  obtain 
the  flraction  of  activity  per  unit  area,  using  ffj^Aj^(r)  to  express  the  distri¬ 
bution  of  activity  in  the  cloud.  Let  us  first  consider  calculating  the 
fraction  for  a  point  on  the  positive  x-axis  at  a  distance,  x  >  R.  If  we 
perform  the  integration  in  the  r-S  plane  indicated  in  Eq.  {'}),  we  wish  to 
Integrate  over  the  area  indicated  in  Fig.  7-  Substantial  simplification  in 
the  final  expression  for  the  fraction  seems  possible  if  we  consider  a 
transformation  from  the  H— r  plane  into  the  t-r  plane,  where  t  is  the  time 
from  stabilization  of  the  cloud  (see  Fig.  8). 

The  relation  between  T,  H,  and  r  is  given  by  Bq.  (V),  where  for  the 
moment  we  consider  q(r)  and  e(r)  to  be  constants.  Thus  our  transformations 
are 

T  -  a 

r 

and  (l6) 

r  ■  r 

The  region  over  which  we  desire  to  Integrate  becomes  that  shown  in  Fig.  7. 

Our  element  of  eurea  is  transformed  by 


d  (H.rl  .  e! 

d(T,r)  S 


(18) 


Fig.  7  —  Schematic  diogram  of  the  region 
of  integration  in  the  r-H  plone 
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'ig.  8  —  Schematic  diagram  of  the  region 
of  integration  in  the  r-t  plane 
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The  integration  of  Eq.  (5)  becomes 


Tg  r2(T) 


vhere 


F(x,y)  =  /  I  A  (r)  drdT 

^  \  r^(T)  ^  P 


(19) 


rwr 

”v  T  -  a 


(20) 


and 


X  +  R 
_ e 

V 

X  -  R 


(a) 


Inserting  the  log-normal  distribution  in  Eq.  (l9)  for  Aj^(r),  we 
obtain 


F(x,y)  = 


PH 

n  _ i 

2  T-a 


// 


Making  the  substitution 

r,  =  ULIHL. 


-i/2(ifLpi) 


^^drdT 


(22) 


1  T-a 


0 


(25) 
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Eq.  (22)  becomes 


-1/2  +  2(qa-Hi) 

e  dq  dT 


where 

Hf  -  a  ^ 

0 

and 

n."“J 

T  -  a  ^ 

0 

(24) 


(25) 


By  coD^letlog  the  sqxiare  In  the  exponent,  Bq.  (23)  becomes 

F(x,y) 


2l«2o^  ^ 

1 - !  I  i. 

^  Tj^  n  .  /  2  TT 


y..  f  f  -l/2(n-2a)‘ 


dT}  dT 


(26) 


Ihe  cumulative  normal  distribution  function,  0(x),  Is  defined  as 
X 


2 

1  r 

9((x)  -  -  J  e  dq 

Then  Eq.  (26)  may  be  written  simply  as 


(27) 


F(x.y) 


\  e 


2u+2a 


/ 


p(m-2a)  -  0(n-2a)3  dT 


(28) 


To  perform  the  integration  with  respect  to  T  we  resort  to  integration 
by  parts.  Since  the  integrations  of  j^(m-2a)  and  j^(fr-2a)  are  functionally 
the  sane,  we  will  perform  only  the  integration  of  0(m-2a)  and  write  down 
the  integral  of  0(n-2a)  by  inspectico.  To  simplify  notation,  let 
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“  =  [-  ln(T  -  a)  +  CjJ 

(29) 

“  “  is  [-  +  s] 

(50) 

where 

Cj  »  In  p  Hj  -  2ti 

amd 

Cg  =  In  p  Hg  -  2ti 

After  making  the  change  of  variables  Indicated  in  Eq.  (29)  ve  obtain 


ni„ 


I  -  /  ji(m-2o)dT  =  (-  2o)  e  ^  /  0(m-2o)“' 


•2om 


dm 


(51) 


>diere 


and 


mi  =  i  in  (Ti  -  a)  +  CgJ 


(52) 


Integrating  by  parts,  Bq.  (51)  becomes 


ran 


,  n  'T  r  ^(nt-2a)e  1  f  -l/2{nh-2of  -2om  ,  1 

I  -  -  -  T=25T2  ;  «  ^55) 

“l  “l 

Performing  the  integration  and  simplifying  Eq.  (35),  ve  have 


I  a  e 


I  ^(mg-; 


-2cmi2  -2cnn^  -2a2 

2o)e  -  0(nL— 2a)e  -o 


|9^(m2)-^(m^)j  I 


(5U) 
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From  Eq,.  (34)  we  can  immediately  write  F(x,y)  as 


F(x,y) 


^(ing— 2a)e  —  0(mj^-2a)e 


-2an_  -2cm, 

0(n2-2a)e  -  ^in^-2o)e 

-e~^°  j^Cng)  -  0(n3^)J 


Inserting  the  values  of  C^,  C^,  n^,  and  n^  in  Bq.  (35)  yields 


2u+2a 


F(x,y)  -  \  ^ 


(Tg-a)  5l(in2-2a)  -  (T^^-a)  0(in^-2a) 
-(Tg-a)  ^(n^-Za)  +  (T^-a)  0(n^-2a) 


To  facllitatf  computation  of  Bq.  (56)  we  have  computed  the  function  ^  (T) 


which  Is  defined  as 


^  (T) 


(T-a)  e 


2u+2cr 


Jp(iB-2a)  -  J^(n-2cT;^-  Hj0(m)  + 


We  may  now  write  F(x,y)  as 


P(x,y)  -  S4  (tg)  -  \f/  (tj^)J 


(t)  has  been  conqputed  on  the  JOHHHIAC  using  for  0(x)  the  approximation 
derived  by  a  suitable  transformation  of  the  function  given  by  Hastings. 


A  local  RAND  digital  computer. 
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Throughout  the  integration  a(r)  and  g(r)  have  been  treated  as  constants. 

It  was  possible  to  plot  the  time  of  fall  vs  a(r)  and  p(r)  with  isopleths 
of  constant  height. In  computing  (t),  a(r),  and  p(r)  should  vary  in 
a  continuous  fashion  from  the  values  on  the  isopleth  for  to  those  on 
the  isopleth  for  H,j,  for  the  particular  T  being  considered.  Since  we 
held  a(r)  and  p(r)  constent  in  ccmputing  \j/  (t),  the  values  of  a(r)  and 
p(r)  corresponded  in  one  case  to  those  on  the  isopleth  for  and  in  the 
other  to  those  on  the  isopleth  for  H^.  The  resulting  curves  for  ^  (T), 
using  these  two  sets  of  a(r)  emd  p(r),  were  the  same  except  for  a  slight 
difference  at  very  small  times.  This  conclusion  partially  Justifies  the 
simplification  of  holding  a(r)  and  s(r)  constant  during  the  integration 
to  obtain  ^  (T). 

Since 

Tg  -  T  =  ^  and  -g-  ■  -  -  —  (59) 

2  1^  2  V 

F(x,y)  may  be  approximated  by 


F(x,y)  =  \  y  '^'(7) 


(40) 


The  derivative  of  ^  (T)  can  be  computed  directly  by  considering  the 


equation 


\l/  (T)  +  C 


f  / 


^/^  ySTa  r 


^  dr  dn 


(41) 


T-a 


where  C  is  a  constant.  Differentiating  Eq.  (4l)  with  respect  to  T  we 


obtain 
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T-a 


f 

w. 


T-<X 


£  ^ 
p 


(U2) 


Substituting  Bq.  (23),  Bq,.  (!(-2)  becmes 


m 


'(T) 


2U 


If:!  / 


.l/2(Ti-2a)' 


dn 


(»^5) 


or 


'(T) 


^2iA+2a 


m 


has  also  been  computed  on  the  JOaNHIAC  for  heights  appropriate  to  yields 
from  l/k  to  l6  MT.  These  curves  are  shown  in  Figs.  9—12. 

Once  the  fraction  has  been  calculated  for  a  point  on  the  x— axis,  it  is 
a  simple  matter  to  obtain  the  value  for  a  point  off  the  axis.  For  a  point, 
(x,y),  the  effective  radius  of  the  cloud  becomes  Rgf(y),  where 


Cis) 


Thus  Eq.  (UU)  becomes 


P(x,y)  -  \  ~  i“(y) 


(46) 


per  unit  time  for  a  16  MT  device 
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Wl- 


S>kf 


Ul 


D.  f(y)  1  -  (5/15)^  -  -9*0 

E.  F(x,y)  =  5*55  X  10~^  x  3  x  .9U3  x  1-57  =  2.1*6  x  10~*^  per  (n  mi)”^ 

To  compute  the  fraction  at  points  with  an  x-coordinate  less  than  we 
proceed  in  a  manner  similar  to  that  above  but  we  use  Bq,.  (1*7)  in  place  of 
Eq.  (46). 
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VI.  DOSE  RATE  AND  IMTBSRATED  DOSE 

Given  the  fraction  down  at  a  point  and  its  mean  arrival  time  it  ia 
poeeible  to  estimate  the  dose  rate  at  any  time  and  the  integrated  dose  to 
any  time.  Assme  that  3  x  10^  gamma  Me  are  produced  per  MT  of  fission 
yield  at  one  hour  and  that  the  hypothetical  dose  rate  at  one  hour  is 
equal  to  3  r/h  per  gamma  Me  per  sq  n  mi.  The  total  number  of  Me  at 
one  hour  produced  by  a  weapon  of  W  MT,  where  f,  the  fraction  of  the 
yield  due  to  fission,  is  3Wf  x  10^.  Hence  3Wf  x  10^  is  the  nmber 
of  Me  per  sq  n  ml  and  the  hypothetical  one  hour  dose  rate  is  given 
by 


-  9WfF  X  10^  r/h  (40) 

•  —1.2 

The  dose  rate  at  some  later  time  is  R^^t  provided  t  is  larger  than  the 
time  at  which  the  fallout  at  this  point  is  coii5)lete.  To  coiiq).ite  the  inte¬ 


grated  dose  at  a  point  (x,y)  we  make  the  assumption  that  ail  of  the  activity 

arrives  at  this  point  at  a  time  equal  to  —  .  Hence  the  Integrated  dose, 

V 

R(T),  at  a  time,  T,  is  given  by 


■  -I  / 


For  convenience  of  conq)utatlon 


i 


RM-2U60 

uu 


5 


(t) 


-0.2 


-  T 


.-0.2 


(51) 


is  plotted  on  Fig.  13  for  the  times  of  usual  interest. 

It  should  be  pointed  cut  that  this  integration  replaces  the  tr’A-^  time 
of  arrival,  vhlch  is  spread  over  an  increment  in  time,  with  an  arrival  time 
appropriate  to  the  arrival  of  the  particle  ftrom  the  center  of  the  base  of 
the  cloud.  Thus  the  dose  rate  will  build  up  during  the  period  of  arrival 
and  fall  off  after  the  rate  of  arrival  becomes  smaller  than  the  rate  of 
decay  as  shown  in  Fig.  lU.  The  arrival  time  used  in  these  calculations 
is  somewhere  near  the  midpoint  in  time  between  the  time  of  first  arrival,  t^, 
and  the  time  of  cessation,  tg.  Cie  true  integrated  dose  is  the  integral 
under  the  solid  curve;  the  integrated  dose  computed  by  the  above  method  is 
the  hatched  area.  The  differences  are  minor. 
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Fig.  13  —  Ratio  of  integrated  dose  to  Ihr  dose  rote  os  o  function  of  near 
time  of  arrival  for  several  final  times 


Dose  rate 
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l(  tj 

V 

Time  ^ 


Fig.  14 —  Estimated 
with  true 


integrated  dose  compared 
integrated  dose 
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VII .  CONCLUSIONS 

Hie  tvo  solutions  of  the  mathematical  fallout  model  provide  a  basis 
for  the  construction  of  simple  and  rapid  computation  of  fedJ.out  contours. 
A  study  now  under  way  Indicates  that  a  slide  rule  may  be  constructed  to 
solve  Bq.*  (Uo).  Another  study  Indicates  that  wind  shear  and  the  vertical 
distribution  with  height  can  be  taken  Into  account  by  applying  the  model 
to  layers  of  the  cloud. 

The  computational  procedures  outlined  here  provide  a  fairly  precise 
solution  to  the  proposed  mathematical  model.  Conputatlonal  procedures  do 
not  Introduce  errors. of  more  than  3  to  10  per  cent.  Hie  degree  to  which 
the  model  reproduces  the  natural  occurrence  may  be  wrong  by  factors  of 
2  to  U.  Improvement  In  the  model  and  the  values  of  the  parameters  In  the 
model  must  wait  for  more  and  better  observations  of  the  phenomena  of 


fallout 


